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ABSTRACT 

We present an analytic description of turbulent, magnetohydrodynamic (MHD) disk accretion around 
black holes that specifically addresses the relationship between radial and vertical, mean-field transport 
of mass, momentum and energy, thereby complementing and extending numerical simulations. The 
azimuthal-vertical component of the magnetic stress is fundamental to an understanding of disk- 
corona-outflow coupling: when it is important for driving the angular momentum transport and mass 
accretion in the disk, it also has an important influence on the disk-corona-outflow energy budget. 
The Poynting flux derived from the product of this term with the Keplerian velocity also dominates 
the Poynting flux into the corona. The ratio of the coronal Alfven velocity to the Keplerian velocity is 
an important parameter in disk-corona-outflow physics. If this parameter is greater than unity then 
energetically significant winds and Poynting flux into the corona occur. However, significant effects 
could also occur when this parameter is much less than unity. A limiting solution describing the case 
of angular momentum transport solely by the vertical-azimuthal stress has the property that all of 
the accretion power is channeled into a wind, some of which would be dissipated in the corona. More 
realistic solutions in which there is both radial and vertical transport of angular momentum would 
have different fractions of the accretion power emitted by the disk and corona respectively. These 
results have important implications for existing accretion disk theory and for our interpretation of 
high-energy emission and nuclear outfiows from the central engines of Active Galactic Nuclei and 
Galactic Black Hole Candidates. 

Subject headings: galaxies: active, quasars — hydrodynamics — magnetic fields 



1. INTRODUCTION 

Since the foundations we re laid for a standard theory of disk accretion l|Novikov &: Thorne| 1973|:|Prin^le fc Rees|1972|: 
iShakura fc Sunvae vI 119731) . two fundamental problems were immediatelv recognized "i^Bi snoTC^^^Cogg^^^Bl^m^ovl 
^77; Liang & Pricd llQTTt IPaczvnski 1978): (1) The efficient transport of angular momentum to large radii cannot 
be attributed to conventional kinematic viscosity, and (2) The observed high-energy spectra and luminosities and the 
ubiquity of outfiow phenomena from accreting black holes implies efficient vertical transport of energy from a relatively 
cool, dense disk to a hot, tenuous and unbound corona. While it has been widely accepted that magnetic fields provide 
the most plausible means of efficiently transporting both angular momentum and energy, the precise nature of this 
transport has remained unclear until recent numerical simulations demonstrated, unambiguous ly, that ac cretion disks 
work because of magnetohydrodynamic (MHD) turbulence (Balbus & Hawlev 1998; Hawlcv. Balbu s, fc Wint ers 1999). 
The turbulence is generated by the magnetorotational instability (MRI - Balbus fc Hawlev. 199ll and references 
therein), which is driven by the free energy available from the differential rotation of the bulk flow. Notwithstanding 
these groundbreaking results, however, the nature of vertical energy transport from an accretion disk to a corona 
and/or outflow still remains an outstanding and contentious issue. 

In the context of Active Galact ic Nuclei (AGN), much theoretical e f fort has recently focussed on the 
physics of accretion disk coronae llDi Matteo. Blackman fc FabiaiJ Il997al: iDi Matteo. Celotti fc Fabiai] Il997hl: 
ILiu. Mineshis-e fc ShibatE]l2nol iMerioni fc Fabianll2n01albl . I200ir commensurate with the dramatic increase in both 
the quality and quantity of high-energy observational data. However, there has been little improvement in coupled 
disk -corona rnodels since the first phenomen ological descriptions of Haardt fc Maraschi (1991, 1993). Current models 
(e.g. ILiu et al.l 2002: Merloni fc Fabian 200^ simply replace the fraction of accretion power transferred from the disk to 
the corona with a Poynting flux quantity estimated from a mean-field buoyant velocity and an equipartition, mean-fleld 
magnetic energy density. While numerical models (e.g. iMiller fc"St onc 200(3) do indeed show that turbulent fluctu- 
ations in a vertically stratified disk are capable of driving the magneto-gravitational modes of the Parker instability 
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()Parkeilll955j) , whether magnetic buoyancy can supply the corona with sufficient power to explain the observed high- 
energy emission is questi onable. Numerical simulations indicate that rna-gnetic buoyancy is an ineffective saturation 
mechanism for the MRI l|Brandenburg et al.lll995l iMiller fc Stondl2000l iStone et allll996D . while the oretical m odels 
for disk coronae require implausibly ideal buoyancy conditions and limiting accretion conditions fe.g. lMerlonil l20031 . 
Realistically, the growth of the unstable buoyant-Parker modes, which is essentially a wave-fluctuation resonance in- 
teraction, must compete against particle-fluctuation interactions, which correspond to dissipation of the turbulence 
and internal heating of the disk. 

The production and ubiquity of outflows from accretion disks around black h oles also remains a challenging prob- 
lem and it is unclear from most theoretical models (e.g. Blandford fc Pavn^ 1982; Li. Chiuch & Bcaclman 199^ 
lUstvQgova et al.l 120001: iLove lacc. Berk, fc Contopolous 1991; Wardlc & Konigl 1993) whether MHD disk turbulence 
plays a significant role (but sec Heinz fc Bcgelman 200 0). Nevertheless, numerical simulations of turbulent MHD ac- 
cretion disks do in fact show that outflows become important in the innermost r egions of turbulent accretion disks (e.g. 
iHawlev fc Balbusll2002HHawlev. Balbus. fc Sto^l200lHStone fc PringOl200lD . Unfortunately, the numerical models 
are restricted by their approach: non-Keplerian motions are defined a priori as fluctuating quantities, so that vertical, 
mean- field transport is not self-consistently taken into account. Indeed, the outflows that emerge in some of the simu- 
lations (e.g. iHawlev fc Ba lbus 2002') are defined as regions where the nett radial flow is outward, rather than inward. 
Furthermo re, the neglect of vertica l angular momentum transport restricts the radial inflow to substantially subsonic 
speeds (seelBalbus fc Hawl evll99^. Inevitably, the resu l ting accretio n rates in the numerical models are typically very 
low (see e.g. Hawlev fc Bal bual2002t IHawlev et a'l]l20nil IStone fc Pri nglc 2001). Numerical models are also restricted 
by computational limitations: simulations that are global as well as vertically stratified are required to estimate the 
fraction of accretion power that can be vertically transported and this is not only computationally prohibitive, but 
also sensitive to numerical dissipation effects, which are difficult to quantify. 

Thus, the present status of black hole accretion disk theory is that there is currently no formalism which self- 
consistently couples turbulent, MHD disk accretion with a magnetically- dominant corona in a framework that can 
accomodate a range of radial inflow and vertical outflow solutions. 

In this paper, we present the first fully analytic description of a turbulent MHD accretion disk coupled to a corona, 
self-consistently taking into account both vertical and radial mean-field fluxes of mass, momentum and energy. By 
deriving the relevant transport equations from first-principles, our formalism provides a non-phenomenological and 
non-empirical approach to the problem of energy transport to a corona and the associated outflow and the inter- 
relationship of this transport with accretion onto the central black hole. In this treatment, we focus on the effects 
of turbulent magnetic stresses and the mean magnetic flux density is assumed to be zero, for the sake of clarity. 
However, our formalism lends itself naturally to the inclusion of nonzero mean magnetic fields, which we intend to 
explore separately in a subsequent paper. In § 2, we present the relevant conservation equations for a resistive, viscous, 
and optically-thick MHD gas. In § 3, we statistically average these equations and derive the corresponding mean-field 
equations for a turbulent MHD gas. In § 4, we apply these mean-field equations to the dynamics of a geometrically-thin 
accretion disk that is stationary and axisymmetric in the mean, and we expressly examine the implications of vertical 
mean-field transport for the conservation of mass and momentum. In § 5, we utilize the results of § 4 to analyze the 
total disk energy budget. We conclude with a discussion of the main results in § 6. 

2. FUNDAMENTAL EQUATIONS 

In this section, we summarize the standard equations for a resistive, viscous, and radiative MHD gas in a dynamically 
important gravitational field. We consider only a nonrelativistic, optically-thick gas in which the radiative diffusion 
approximation holds, so that the radiation pressure reduces to a scalar. The independent variables are: mass density, 
p; fluid velocity, w^; gas plus radiation pressure, p; gas plus radiation energy density, u; gravitational potential, 0g; 
radiative flux, Fi] external heat fiux, Qi] magnetic field, Bi] electric field, Ei] current density, J^; and viscous stress 
tensor, tjj^ 

The fundamental equations are: 



1. Mass continuity. 

1 + ^ = . (2-1) 

Ot OXi 

2. Momentum conservation. 

dv dv, djpv,) , d dcj,^ dp d /B,B, dti 



3. The induction equation for magnetic fiux conservation, 

dB^ _d_ 



OB d 

-KT + ^vkekim-^ — {BiVjn) = ??V^Bj , ( 2-3) 



Note that this is the viscous stress tensor describing microscopic processes. It does not represent the so-called turbulent viscosity. 
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which is derived from Maxwell's equations, 



and a scalar conductivity law, 



Ji^cr (^E, + ^e^jkVjBk^ , ( 2-5) 

„2 , 



with a the conductivity and rj = c /Ana the resistivity. 
4. Internal energy. 



du 3 j'^ 

-TTT + ^(""^i) = -pvi,i - Fi^i - Qi^i H h tlj-Wij . ( 2-6) 



5. Electromagnetic energy. 

We define the magnetic energy density, Ub, the magnetic stress tensor, tf^, and the fluid shear tensor, Sij, by: 

(2-7) 
( 2-8) 





~ Svr 






Bj 
Air 


Sii — 
^ Btt 


Sij 







S^JVkM ] ■ ( 2-9) 



Taking the scalar product of the induction equation, (|2-3|l . with Bi gives 

Oub d 1 1] 



+-Qj^{uuVj) + ^UBVkM^tfjS,j + ^B,V^B, , (2-10) 



This equation describes the volume rate of change of magnetic energy due to advection, expansion or compression, 
shearing and diffusion of field lines. We write this equation in this particular form so as to emphasize the shear 
term, which is the source term for magnetic field amplification in a fluid with shearing motions, such as an 
accretion disk. Equation (|2-l()|l also shows the well-known result that in the absence of shear or in the case of 
an isotropic magnetic field the magnetic energy density evolves similarly to a gas with adiabatic index 7 = 4/3 
(i.e. pressure = 1/3 x energy density). The diffusion term in (|2-10|) satisfies 

-^B.v^B, = ^rir^ ' ( 2-11) 

47r a oxiOXj 

which describes the effects of Joule heating and the diffusion of field lines. Equation (|2-10|) can then be expressed 
as 

5mb d 



dt dxi 
Note that the Poynting flux is 



UbVi + 1]- 



1 J2 

= tfjStj - ■:^UsVk,k • ( 2-12) 



a 



c dt^- 
Si = —eijkEjBk = u^Vi - tfjVj + T]-^ , ( 2-13) 



Si^i — —JiEi — —tijk — JjBk — ( 2-14) 



Equation H2-12|l can then be expressed in the more familiar form 

Oub „ , „ Vi , „ J"^ 

dt - c ' G 

describing conservation of electromagnetic energy. 

All of the above forms of the equations describing the evolution of the magnetic energy in a fluid are useful. 
Nevertheless, it is worth emphasizing that when an induced current is involved in JiEi (as it is here through the 
expression for the current), the Poynting flux does not correctly represent the flux of electromagnetic energy. 
That is, the Poynting flux is only a good representation of the flux of electromagnetic energy in vacuo. This 
point that is not widely appreciate d although it has been known for some time in the plasma astrophysics of 
emission processes (see for example. iMelrose &: McPhedranI l)1991|) V This point assumes some importance when 
we consider the disk energy budget in § 
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6. Total energy. 

We define the specific entlialpy of gas plus radiation by 



h = ^ ( 2-15) 

P 



and we also define the total energy and corresponding energy flux by 

1 
2 



"tot = ^pw^ + P(/>G + u + Ub ( 2-16) 



Ff = Qpi;' + + phj v, + F,+ Q, + 5, - t^^jV, . ( 2-17) 

The conservation equation for the total energy is obtained by taking the scalar product of the momentum 
equation, H2-2() . with Vi and then utilizing both the internal energy equation, l|2-6l) . and the electromagnetic 
energy equation H2-14|l to obtain 

%^ + f^=0 . (2-18) 
ot oxi 

Equation (|2-17|l is the more usual form for the energy flux that explicitly incorporates the Poynting flux. Nev- 
ertheless, a more useful form for our purposes is derived by substituting the expression (|2-13|) for the Poynting 
flux, giving 

/ 1 \ dt^- 

Ff = + + p/j + j V, - tf^vj +F, + Q, - t^vj + . ( 2-19) 

In this expression, the Poynting flux is mainly replaced by an advection term u^Vi and a term t^jVj representing 
the rate of work by magnetic stresses on the flow. 

3. STATISTICAL AVERAGING - THE MEAN-FIELD EQUATIONS 
3.1. Definitions of mean variables 

In t he follow ing, we make extensi ve use of a statistical averaging approach (e.g. lBicknelllll984. 119861 lBradshawlll976t 
lFavrel ll969: K rause fc Radlerlll980tl in which the MHD equations are ensemble-averaged. When the flow is steady in 
the mean, we can think of the averaging process as involving an average over a time-scale large compared to the time 
scales of instabilities and turbulent fluctuations. Statistical averaging is not commonly used in astrophysics ( at least 
not explicitly). However, it is commonly used in the theory of the physics of turbulent fluids and in application areas 
such as geophysical fluid dynamics. The approach is capable of providing valuable insights. For example, in the case 
of accretion disks we gain valuable insights into the energy flow within the system. 

To describe the mean flow dynamics, we adopt the prescription whereby a quantity, X, is expressed in terms of 
mean and fluctuating components and {X) denotes its ensemble average value. We distinguish intensive quantities, 
such as velocity, from extensive quantities, such as density, pressure and magnetic flux density, by using mass-averaged 
values, X, and mean values, X, respectively. This has the effect of conserving the mass associated with the mean 
flow, of conserving the mean magnetic flux and of restricting the number of terms that are produced by the averaging 
procedure. Thus, the statistical averages of the density, velocity, pressure and magnetic field are represented by: 

P - P + P' ip') - 

p = p + p' {p') = ^ '^-^> 

R = B, + B' IB') = 



and the mean fluid shear is 



\ + - '^^^jVk,k ] ■ ( 3-2) 



In this paper, we restrict ourselves to the simplest case in which the magnetic field is dominated by its fluctuating com- 
ponent, so that {B) = B <C (B''^)^!'^ . The more complex case in which a systematic component of the magnetic field is 
present clearly has important implications for disk accretion solutions, particularly when the field strength is sufficiently 
high to quench the MRI and transport angular momentum via larg e-scale m agnetic torq ues such as those invoked in 
models for MHD-driven outflows and Poynting flux jets (e.g . Blandfor^^_^P|3,vne 1982; L i. Chiueh &: Begelmanl[l992l: 
lUstvQgova et al.ir200ft iLovelace. Ber k, fc Cont oDoloudl ll991t iWardle fc Konig l 1993). However, it is unclear whether 
large-scale, organized mean fields can be generated from a highly chaotic underlying fl ow and whether they can explain 
collimated outflows associated with accretion disks (see e.g. the discussion in Hein z fc BegelmanI l200(y . The more 
general case of 5 is thus defered for future work. AH the mean-fleld equations derived here can, however, be 
straightforwardly generalized to include terms involving B that are directly analogous to their turbulent counterparts. 

The molecular viscous stress tensor can also be written as tj'^ = i^^ -I- t^-^ , with (i^ij) = 0, since it can be related 
to a coefficient of kinematic shear viscosity, i/, by fl^ = 2ypsij (ignoring bulk viscosity) , where is the shear tensor 
defined in H2-9() . In practice, the mean viscous stress, {tj^) = tjj, usually has a negligible effect on momentum transport 
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(particularly in a high Reynolds number turbulent flow) and is thus usually ignored. However, the fluctuating part 
of the viscous stress plays an important role in the dissipation of turbulent energy on the smallest scales at the end 
of a turbulent cascade. Mathematically, this is described by the appearance of the correlation term {tjjvl j) in the 
mean- field internal energy equation, (|3-13|l . derived below. This term represents the mean rate of turbulent viscous 
heating and is dominated by the high wavenumbcr components of the turbulent velocity fluctuations; that is, by the 
dissipative region (in wavenumber space) of the turbulent cascade. Thus, although we henceforth assume (tj^) = 0, 
we retain tjj throughout the averaging procedure in order to consistently include stochastic energy dissipation by 
turbulent viscous stresses. 

3.2. Mass conservation 
Statistical averaging of the continuity equation, (|2-l|l . gives 

dp d{pvj 



dt dxj 

so that the mass defined by the mean density is conserved. 



= , ( 3-3) 



3.3. Momentum conservation 
In the case of zero nett magnetic flux, the magnetic stresses are 

B'B' ^ B'2 



with ensemble mean, 



We also define the Reynolds stresses as 



t"^. ^ ' J ^ (5 -^^ ('3-4^ 



(t-)-^-^-S-!^ (35) 



tfj=-{pvlv'^) (3-6) 
Statistical averaging of the momentum equation, H2-2|) , yields 

This equation shows that momentum transport in a fiuctuating MHD fiuid involves additional, statistically-averaged 
quantities that arise from the turbulent magnetic stresses and the Reynolds stresses. Note that the Reynolds stress is 
conventionally defined as a mean quantity, whereas, for convenience, we have defined tfj as an unaveraged quantity in 
order to simplify triple correlation terms that appear in the energy equations derived in § | 3. 51 b elow. Note al so that 
these definitions of the turbulent stresses are opposite in sign to those used, for example, bv lBalbus fc HawlevI lfl998i) . 
Our usage is consistent with the conventional description of this term as a stress. 

Often in astrophysical contexts, the terms tfj and tfj are referred to as viscous and resistive terms since they represent 
additional forms of momentum and electromagnetic transport. In this paper, we consistently refer to these terms as 
turbulent stresses and reserve the term viscosity and resistivity for real molecular effects. 

3.4. The induction equation 

Although we assume that the mean magnetic flux density is zero, for completeness, and to make one physical point, 
we present here the mean-fleld induction equation obtained by statistically averaging H2-3|l : 

dB- 3 - - 

-g^ + eyfeCfc;™^ {Biv„, + {B[v'J) = yyV^B, ( 3-8) 

This equation forms part of the foundation of classical mean-fleld electrodynamics (e.g. lKrause fc R,adledlT98?1) . which 
is usually developed under the assumption of an incompressible flow. The contribution from the term involving Biv„i 
results from compressibility, while the term involving (B ' yl„) g ives rise to dynamo ampliflcation resulting from mean 
helicity. As exceptionally clarified bv lBalbus fc HawlevI l)1998l see Sec. VI-B), when a weak seed field B is present 
in a differentially rotating ffuid, the fluctuating flelds v' and B' are self-consistently generated by the MRI and thus, 
cannot be prescribed independently of B, which is the underlying assumption of mean-fleld dynamo theory. We do 
not, therefore, make any assumption in our formalism concerning mean-field dynamo amplication. 

3.5. Energy 

There are a number of different energies to account for in a turbulent MHD fiuid: the turbulent magnetic and 
kinetic energies, the internal energy and the total (turbulent plus internal plus mechanical) energy. Here, we derive 
the conservation equations describing the evolution of these quantities. 
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3.5.1. Turbulent magnetic energy 

The mean-field equation for the turbulent magnetic energy density, (ub) — {B''^)/8Tr, is derived by statistically 
averaging the electromagnetic energy equation, (|2-12|l . yielding 

^(ub) + ^ [Mv^ + M] = {t-)S,, + {t-s',^) huu)vk,k iM^k) - Vp§^ - ^ ' ( 3-9) 

where the turbulent magnetic stress tensor, tfj is defined by H3-4|l . The first term on the right hand side of H3-9|l is a 
source term for the production of turbulent magnetic energy as a result of interaction between the turbulent magnetic 
stresses and the mean fluid shear. The second term represents the interaction between the turbulent magnetic stresses 
and the fluctuating components of the mean shear. The third and fourth terms describe changes in the turbulent 
magnetic energy due to compression and expansion in the mean and fluctuating components of the flow, respectively. 
The last two terms on the right hand side of (|3-9|l describe the total work done by resistive forces minus the rate at 
which turbulent magnetic energy is converted into heat via Joule losses. 

3.5.2. Turbulent kinetic energy 

A similar mean-field equation for the turbulent kinetic energy density, (mk) = {■^pv'^), can be derived by taking the 
scalar product of Vi with the unaveraged momentum equation, (|2-2|l . statistically averaging this equation and then 
subtracting from it the scalar product of Vi with the statistically averaged momentum equation, l|3-7|) . In symbolic 
terms, this is (v • momentum eqn.) — v • (momentum eqn.), which is equivalent to (v' • momentum eqn.) and yields 



<9 , > d 



Uti< 



— [(^iK>«. + M) + (pvl) - {tf^v'^) - {t^jV'^)] 



dt dxi 

= iSS»j-^("K>%fe-(i^4) + (^"B<fe> + (K.)-(*J}<j) ■ (3-10) 

Note the similarity to the turbulent magnetic energy equation, (|3-9|l . and also to the internal energy equation for an 
adiabatic 7 = 5/3 gas when the triple correlations are neglected. In particular, there is an analogous source term, 
tfjSij, describing the rate at which shear in the mean flow does work on the Reynolds stresses. There is also a sink 
term, {t'^jV^ j), describing viscous dissipation of turbulent kinetic energy into heat. Note also that the turbulent kinetic 
energy equation is coupled directly to the turbulent magnetic energy equation via the magnetic term on the right hand 
side of (|3-lUf) . which also appears explicitly on the right hand side of (|3-9|1 . 

3.5.3. Total turbulent energy 

The magnetic coupling term in (|3-10|l implies that the turbulent kinetic and magnetic equations can be combined 
into a single turbulent energy equation: 



d . .8 



d{t. 



dxj 



2 1 (J2) 
■ {Uj)~Si] - (gWK + -UB)wfe,fe + (pu-^,) - (ty-w^j) — , ( 3-11) 



where 



(t,,) =tf^ + {tf^) ( 3-12) 

is the combined Reynolds and magnetic stress tensor. 

The last two terms on the right hand side of (|3-ll|l represent the dissipation of turbulent kinetic and magnetic energy 
whereas the flrst term on the right represents the production of turbulent energy through the action of the mean shear 
on the total stresses. The usual approach in many physical applications is to equate the heating rate to the production 
term. However, the presence of transport terms, as well as terms describing the work done by fluid compression and/or 
expansion, implies that the steady-state production and dissipation of turbulent energy are non-local and therefore, 
the rates cannot in general be equated. Thus, it is not generally correct to simply replace the viscous dissipation 
term in the internal energy equation with the production term for turbulent energy. The rate of turbulent viscous 
dissipation in the mean-field internal energy equation must be related to the rate of production of turbulent energy 
self-consistently using H3-ll|l . as is demonstrated in the following section. 

3.6. Internal energy 
Statistically averaging the internal energy equation, (|2-6ll . yields 

8Ti 8 P 

^ + ^ {Ui^ + {UV'^) = -pi^^, ~ {PVU) - ~ (Qm) + (" ) + {t^A,) ■ ( 3-13) 

The terms on the left hand side describe the total rate of change in the gas plus radiation internal energy density 
as a result of intrinsic temporal variations and advective transport by the mean plus turbulent fiow. The terms on 
the right hand side of (|3-13|l describe the work done by compression or expansion in the fiow against the gas and 
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radiation pressure, radiative losses, energy exchange from an external heat source or sink, mean-field Ohmic heating 
and turbulent viscous heating. The source terms determine the rate at which energy is converted into random particle 
energy (some of which is then converted into radiation) and also into bulk kinetic energy. Magnetic energy in particular 
can be converted directly via Joule heating (usually identified with field line reconnection) as well as via work done 
by the turbulent flow against the turbulent viscous stresses 

3.7. Turbulent plus internal energy 

The viscous and Joule dissipation terms appear as sink terms in the turbulent energy equation (|3-f l|l and as source 
terms in the internal energy equation l|3-13|) . Hence, combining these equations eliminates these terms with the 
following result for the total internal energy of a turbulent MHD fluid: 



d_ 
dt 



b)) 



d 

dxi 



(ph 



+ (Wk + Wb) i^ + (ph'v'i) + ((liK + Ub) vl) ~ {{tfj + t^) v'j) + 77 



d(tf 



{U 



1 



( 3-14) 



3.8. Total energy 
Statistically averaging the total energy equation, (|2-18(l . gives 

dutot , d{Fl 



where the average total energy is: 



dt 



"tot = -^PV^ 



dxi 



= 



( 3-15) 



( 3-16) 



and involves both the components that one expects from the mean flow plus the average turbulent kinetic energy, 
(uk)- We use the second form of the energy flux, equation (|2-19() . to form the mean energy flux: 



1-2 

-^pv 



ph 



-{ph'v[) + (mkU-) + (ubw: 



{tij) 
dxi 



( 3-17) 



The numerous terms in this expression for the mean energy flux, (Ff ), are all easily interpreted. The first group of 
five terms represents the energy flux advected by the mean flow and consists of bulk kinetic, gravitational, enthalpy, 
kinetic and magnetic terms. The next group of three terms {{ph'v[) -\- {u^v[) -\- {u^v[)) represent the turbulent fluxes 
of enthalpy, turbulent kinetic energy and turbulent magnetic energy. The next two terms {—{tij)vj — {tijv'^)) represent 
the nett work done by the total Reynolds plus magnetic stresses on the mean and turbulent flow. The next two terms 
{{Fi) + {Qi)) represent the radiative and heat fluxes respectively. 

The last two terms (— {tJjVj) + iv^tfj) j) represent the work done by the viscous and resistive stresses. In most regions 
of a high Reynolds number flow, they are negligible. Nevertheless, in some regions of high spatial gradients, such as a 
shock or a reco nnection region, they could be comparable to the other terms and their inclusion in equation H3-17|l is 
logical (see e.g. iLazarian fc Vishniadll999l for further details on stochastic reconnection). However, when the energy 
equation is integrated over a volume whose bounding surface is well outside regions of dissipation, the contribution of 
these two terms to the resulting surface integral of the energy flux is negligible. Therefore, in many circumstances, we 
can safely ignore these terms in considering the integral form of the energy equation. This point is elaborated in § |21 
in the context of accretion disks. 

The total energy equation describes the nett transfer of energy from one component to another and incorporates the 
work done by the total turbulent stresses through the term {tij)vj as well as the change in binding energy represented 
by the terms ^pv^ + p(t>G- The volume dissipative terms {{t^jSij) and (J^/cr)) disappear because of their equal and 
opposite contributions in the equations describing production of turbulent kinetic and magnetic energy on one hand 
and dissipation of that energy into heat on the other. Other terms in the energy equation represent non-local effects 
such as advection and diffusive transport. 

Note that the advection terms contain contributions from both turbulent kinetic and magnetic energy, not just 
enthalpy. Hence, in applications to advective accretion disk models, for example, the magnetic and turbulent kinetic 
energy should be taken into account especially when the magnetic field is near equipartition or when the turbulent 
velocities are near transonic. 



3.9. Poynting flux 

For future reference (in § \^ , we also note the mean Poynting flux: 



( 3-18) 
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3.10. Mechanical energy of the mean flow 

For completeness, we note that an equation for the mechanical energy can be derived either by taking the scalar 
product of the momentum equation, H3-7() . with the mean velocity Vi, or by subtracting the total internal energy 
equation, H3-14|l . from H3-15|l . giving 



d 


'1 -2 




d 








at 


2^^ + 


P<f>G 











dp d{t. 



^ dxi dxj 



= . ( 3-19) 



4. ACCRETION DISK DYNAMICS 



We now apply the generalized mean-field equations derived in the preceding section to an accretion disk around a 
black hole. In this section we consider the implications of the radial, azimuthal and vertical momentum equations. We 
then apply these results to the energy budget in acccretions disk in § |S1 

In the following, we present the full statistically-averaged equations in a cylindrical (r, (jj, z) coordinate system for 
a fluid which is time-independent and axisymmetric in the mean {{d/dt) — {d/d<j)) = 0). We use Newtonian physics 
throughout, with the gravitational potential 0g — —GM{r^ + z^)^^/^. We integrate the mean-field conservation 
equations vertically over an arbitrary disk scaleheight, h = h{r). Quantities calculated at the disk surface (z = ±h) 
are denoted by a ± superscript, X^, and we assume refiection symmetry about the disk midplane, so that \X^\ — \X^\. 
Midplane values of variables are denoted by Xq. 

Since we include the effects of a disk wind, we do not identify /i as a hydrostatic scaleheight, but as a photospheric 
height, delineating between the disk proper and the transition region leading to a corona, by analogy with the solar 
atmosphere. We assume h is much less than the radius, so that quantities of order h/r and dh/dr are neglected. The 
following comments are useful in the context of the vertical integration process. Let us represent a generic conservation 
law for mass, momentum, energy etc. by 

IdirAr) , OA,, ^ 



dr dz 



4-1) 



Multiplying by 27rr and integrating with respect to the vertical coordinate, z, from —h to h gives: 

y +^ ^ 2^^^^+ - a; ) - y 27rr S dz ( 4-2) 



Taking the derivative with respect to r outside of the integral gives: 

r+h 



d_ 

dr 



n dh 

/ 2TTrArdz-2nr{A++A:;)—+2TTr{A+~A^)^ 2nr S dz 
J-h ar J_h 



(4-3) 



Since reflection symmetry about the midplane holds, 

|A+| = |A-| and \At\ = \A;\ . (4-4) 

Generally, there are two surface terms resulting from the integration over disk height, the second and third terms on 
the left of equation l|4-3|) . Unless \Ar\ >> \Az\, the flrst surface term, which is proportional to dh/dr = 0{h/r) is 
negligible compared to the second, so that the result of the integration over z is: 



d 

— / 2TTrArdz + 2TTr(At - A-) ~ / 27rr5dz 
d»-7-h J-h 



(4-5) 



We believe the geometrically-thin disk approximation ijLvnden-Bell fc Pringlelll974|) provides the most physically 
plausible disk solutions, given that geometrically-thick, advection-dominated disks are founded upon assumptions 
(namely preferential ion heating, negligible ion-electron coupling, and n egligible electron heating overall) that are 
highly i dealized, particularly in the presence of MHD turbulen ce (e.g. iBisnovatvi-Koean fc Lovelacd 119971 120001: 
see also lMerloni fc Fabianll200^ . Furthermore, 'Beeelman' (2002) has recently demonstrated that radiation-pressure- 
dominated solutions need not be concomitant with geometrically-thick disks either. 

Since the MRI is a weak-fleld instability, it drives subsonic turbulence, with {pv'^)<{pCg) <C p^k, where Cs — 
(kT / prrLpY^"^ is the local sound speed and = [GMjr'^Y^'^ is the Keplerian an gular velocity. An i mportant distinction 
between our approach and that adopted in numerical simulation models fsee lBalbus fc Hawlevl lT998'l is that we do 
not identify the mean and fluctuating fluid velocity components with the Keplerian and non-Keplerian components, 
respectively. Instead, we adopt the more formal statistical averaging approach (as outlined in Secl^J of decomposing all 
velocity components into mean and fluctuating parts. Thus, we explicitly distinguish between mean radial, azimuthal 
and vertical motions and their fluctuating counterparts and only the fluctuating components are restricted to subsonic 
speeds. The only restriction we place on the mean fluid velocity components is that they satisfy 3> VriVz and that 
and {i^ are independent of z, simplifying the vertical integration. Thus, we are able to self-consistently treat radial 
and vertical transport of mass, momentum and energy by the mean fluid velocity fields. 
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4.1. Mass transfer 
Vertical integration of the mean- field continuity equation, 

1 d{rpvr) d{pvz) 



dr dz 



= 



gives 



2'KrpVr dz + Anrp v]^ = 



We now introduce the usual definitions for the surface mass density, 

r+h 

E(r) = / pdz 
J-h 

and mass accretion rate. 

We also introduce an analogous mass outflow rate. 



Mwir) = / 47rrp+u+dr = Mwin) - / Anrp^v^dr 

<J r •J T\ 



associated with a mean vertical velocity v'^ at the disk surface, i.e. at the base of a wind. 
Using the above definitions, the vertically integrated continuity equation, (|4-7|l . can be written as 



(4-6) 
(4-7) 

(4-8) 
(4-9) 

( 4-10) 



dr 



M^{r) = Awrp+v^^ 



dr 



M^(r) 



implying that 



Ma(r) -I- Mw(r) = Ma(n) -I- My^in) = constant = M 



( 4-11) 



( 4-12) 



where Ma(ri) -I- Mw(n) is the total mass flux at the innermost stable orbit, n. Equation H4-ll(l implies that under 
steady-state conditions, the radial mass flux decreases towards small r at the same rate as the vertical mass flux 
increases in order to maintain a constant nett mass flux, M, which is equivalent to the nett accretion rate at r — oo. 



4.2. Radial momentum 
The radial component of the mean- field momentum equation, H3-7f) . is 



9 f -~2 

Yr ("^"^ 



d . . 
oz 



dp 
dr 



— {r{trr)) - {t4,4,) 



d{tr 



dz 



( 4-13) 



Several of the terms in this equation are negligible compared to the dominant term, pv'^/r. To determine which terms 
can be neglected, we first use the mean-field continuity equation, 1)4-6(1 . to simplify the d/dr and d/dz terms on the 
left hand side. Then, evaluating the gravitational term on the right hand side gives 



l_dvl 
2^ dr 



GMp dp 1 

J.2 Qj, J, 



d{tr 



dz 



( 4-14) 



The radial gradient terms can be dropped since they are negligible compared to pv'^/r. However, the vertical gradient 
in the turbulent stress, d{t^z) / d'^i could be important in a geometrically thin disk, so this term is retained. Integrating 
the remaining terms in (|4-14|l over z gives 



GMT. 



which reduces to 



2r 



2{tr 



( 4-15) 



( 4-16) 



Thus, Keplerian rotation prevails in regions where the turbulent rz stresses on the disk surface are < h^^jr times 
smaller than pw^, where /lav is the density scale height defined below (see equation H4-18f) 'l. Since the turbulent 
stresses saturate at quasi-thermal levels, they are unlikely to modify the Keplerian profile. Interestingly, models for 
flux emergence from the solar photosphere indicate that the Parker instability can enhance the poloidal flux sufhcently 
to modify the background shear in the rotation velocity, thus giving rise to an effective buoyant-shear instability (e.g. 
IChne. Brummell fc CattalIecll2nol . In numerical simulations for accr etion disks , however, this does not appear to be 
the case, with Keplerian rotation profiles emerging in all models fe.g. lHawiev"fc Balbus..20021 . 
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4.3. Vertical momentum 
Vertical momentum balance in the disk is obtained from the z-component of eqn. (|3-7|l : 

We define a density averaged disk height and a density averaged vertical velocity by: 

2 rh 2 ,h 



hav = / pz dz , -D^^av = — pvz dz . ( 4-18) 

^ Jo ^ JQ 

Integration of equation (|4-17l) over z then yields: 

(-Ma^;.,av) + ^ 2r^Po -P + - ' ( 4-19) 

where the rz stress term has been dropped since its contribution is ~ ha,v/r that due to the vertical pressure gradients. 

In order to determine the relative importance of the remaining terms, note that the first term on the left hand side 
of H4-19(l can be written as 

where equations (|4-9|I - H4-11|I . involving the mass fluxes, have been used. The first term on the right hand side of this 
equation is Vz,a.v/vt ^ 1 times smaller than the pw^^ term in H4-19I) and is thus negligible. The relative importance 
of the second term on the right hand side of the above equation depends on the radial gradient of Vz ; if we reasonably 
suppose that dvz/dr ^ Vz/r, then this term is smaller than the pw^^ term in the vertical balance equation, H4-19|l . 
by a factor ^ h/r. We thus neglect altogether the contribution from the first term on the left hand side of H4-19|l to 
the overall vertical momentum transport in the disk. Thus, vertical momentum balance in the disk implies a mean 
vertical outflow from the disk surface given by 



The following well-known order of magnitude estimate for the hydrostatic equilibrium disk height is useful and we 
repeat it here for the sake of completeness. Equation (|4-21|l with -u^ = and E ~ 2po^av (po — central disk density), 



•-av 



r 



( 4-22) 



where cq = \/po7po is the generalized sound speed corresponding to the total central disk pressure Pq°* =po + {ub)o- 
If the departures from equilibrium values of pressure, magnetic field etc. are by no more than a factor of order unity, 
then the velocity at the base of the wind, u+ implied by l|4-21() is of order cq. 

4.4. Conditions for a disk wind 

In the context of the specific accretion disk application being considered here, a non-advective, zero nett magnetic 
flux disk, with turbulence driven by the MRI instability, the thermal pressure cannot drive an energetically-significant, 
large-scale wind, not even for a coronal temerature ~ 100 keV. For example, by examining the Bernoulli equation 
(without magnetic field), it is straightforward to show that the thermal temperature required to drive an outflow in 
the vicinity of a black hole is kT > ^(7— l)GMiimp/r = \v1^^ ~ 10^(rg/r) keV, where 7 is the adiabatic index, Vesc is 
the escape speed, and r^ = GMj (? ~ 1.5 x 10^^M7 cm is the gravitational radius of a black hol e of mass M = 10 ^ M^Mp,. 
If the disk luminosity is super-Eddington then a radiation pressure driven outflow is feasible (iBegelman "2002) . On the 
other hand if the radiation field is sufficiently diffusive to smear out vertical gradients, then radiation pressure may be 
unimportant (^laes 2002:). 

The most probable driving source for a large-scale wind is the magnetic field and one may derive approximate 
conditions for the production of a wind as follows. Let be the Alfven speed, with mean value defined by 

^. = K>^/^==(^)''' (4-23) 
and consider the vertical component of the energy flux at the disk surface without turbulent diffusion terms: 

{F-)+ ^ \^-v+' + ^o + h++ j P+vt - V. . ( 4-24) 

The azimuthal magnetic field component, B'^ should dominate so that we put — — ~ (^^a)- There are negative 
(J^v^ + 'pG — ~5^T^) a-iid positive {v'j)^) contributions to the first bracket of terms in the energy flux. The last 
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term (the centrifugal term) makes a positive contribution to the energy flux if the mag netic stress, (in)-'^ {B'^B'^)+ < 0. 
In a steady state, the energy flux, integrated over the extent of the wind, is conserved and in order that the wind 
escape to infinity, the integrated energy flux should be positive. Neglecting, in the first instance, the centrifugal term 
and disregarding the specific enthalpy, the energy fiux is positive if 



( 4-25) 



that is, the Alfven speed in the corona should be comparable to, or exceed, the local Keplerian speed. 

Now let us suppose that the magnetic stress term is systematically negative, and that owing to the factor of the 
Keplerian velocity, the term —{Att)^^ {B'^B'^)vk, is significantly greater than Then, the condition for the 

energy flux to be positive is: 



4tt 



> U< + P^K ( 4-26) 



and since 0g = ^''^k using equation H4-ll() for dM^/ dr, this condition becomes 

dM, 



dr 



<87rr ^ (4-27) 
47ruK 



that is, the existence of a wind, in this case implies an upper limit on the wind mass-loss rate. 

A more detailed investigation of the conditions for the initiation of disk winds and the corresponding mass-loss limits 
is deferred to future work. We note, however, the results of Meier ct al. (1997) and Meier (1999), pertaining to the 
case of a nett poloidal field, showing a change in character of disk winds, from loosely coUimated to jet-like when 
the Alfven speed increases from below to above that parameter. Their results and the above order of magnitude 
estimate suggest that the condition (|4-25|l may be an interesting critical value. Note, however, that equation H4-27|l is 
potentially a weaker condition on the magnetic field, and can be expressed as: 



{vlY 




( 4-28) 



Given that the vertical wind velocity (at the base of the wind) is likely to be no greater than the sound speed and 
since hgy/r << 1 for a non-advective disk, the right hand side of this inequality is much less than unity. Thus, even 
allowing for the fact that the vertical component of the magnetic field may be significantly less than the magnitude 
of the magnetic field, the coronal Alfven speed implied by equation H4-28|) should be considerably less than the local 
Keplerian speed. Notwithstanding these estimates and the work by Meier and colleagues, it will be of interest to assess 
the validity of these estimates numerically especially in the case of a zero nett magnetic flux. 

Satisfaction of the above conditions (but especially H4-28|l ) is possible in principle, even if the magnetic field in the 
interior of the disk is weak enough for the MRI instability to apply. If magnetic field is driven into the corona by 
buoyancy and the magnetic field does not decrease as rapidly as the square root of the de nsity then the Alfven speed 
in the corona could become quite high. This is, in fact, a feature of the iMiller fc Ston e (2000) simulations which 
therefore support the feasibility of a disk wind, in the case of zero nett magnetic flux, even though a strong wind was 
not produced in their work. 

4.5. Angular momentum 

The azimuthal component of eqn. is: 

{r-pvri,) + I ipi.i.) = [rHu,)) + ^ . ( 4-29) 
Let us define the vertically integrated r0 stress by 

r + h 

Tr^= / U^dz (4-30) 



J-h 

Integrating equation (|4-29|l over z and applying mass continuity H4-ll|l gives 

d 
dr 



Ahv^r + 2nr^Tr^'\ = ^'1'^'^ - 47rr2(t^,)+ ( 4-31) 



The terms on the the left hand side of this equation describe radial transport of angular momentum due to the mean 
radial infiow and turbulent MHD stresses, while the terms on the right hand side describe vertical angular momentum 
transport due to mass loss in a wind and turbulent stresses on the disk surface. Although we do not explicitly consider 
the effects of a nonzero mean magnetic field here, we note that equation (|4-31|l can be straightforwardly generalized 
to include electromagnetic torque terms analogous to the turbulent MHD stress terms involving T^^ and (i^z). 
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Radially integrating H4-31() gives 



dr 



dr 



and using the integrated continuity equation, Ma + Af^ = M , this can be conveniently expressed in the form: 



where 



C(r) = 1 



dr 



W0(r)r 



( 4-32) 



( 4-33) 



( 4-34) 



Note that we do not assume the r(p stress is zero at r = rj. While the assumption Tr^{ri) ~ is often made in 
conventional disk models, its justification rests on the premise (in a-disk models) that the stresses scale as the disk 
gas pressure, which falls off dramatically as the matter passes accross the marginally stable orbit. When the stresses 
are tu rbulent, on the other hand, this assumption is no longer valid. Indeed, numerical simulations l)Hawlev fc KroliS 
120021) show that turbulent stresses persist at and inside ri. As can be inferred from H4-33|l . a nonzero rtj) stress at 
Ti requires an increased torque at large r to transport the same amount of angular momentum over the same radial 
distance. The fact that Tr^{ri) does not generally vanish also has important implications for energy balance in a 
turbulent disk, since this stress does work on the rotating disk and the resultant turbulent energy can be dissipated 
at a comparable rate, thereby providing an additional source of internal heating at rj. This is demonstrated explicitly 
in §0 below. 

An important property of MRI-driven turbulence is that the rcj) stress has the same sign as dn/dr (IBalbiis fc HawlevI 
Il998j) and thus, even in the absence of a nett vertical angular momentum flux, the Tr^ stress alone can facilitate 
the outward transport of angular momentum required for accretion to proceed. Whether accretion proceeds at an 
interesting rate, however, is another issue deserving separate attention. Consider a Keplerian disk in which the mass 
accretion rate varies as Mf^jr) jx rP, where < p < 1 is the mass loss index used, for example, bvlBlandford fc Bee relmaiil 
((19 9 Q.") (see also .Becker ct aLl l2001l) . Substituting into the integrated angular momentum equation H4-32|l yields 



2p+l 



1 - 



-{p+- 



Gr4,{r) 



1 - 



Gr<p{n) 
Gr4>{r) 



G^z{n) 



G^.{r) 



( 4-35) 



where Gr4,{r) — 27rr^(— T^^) and G^z{r) = — 47rr^dr(t02)+ are the torques associated with the turbulent MHD 
stresses. This implies that if Gr^, G^z cx r^'^^Z'', then 



M^v^r ~{Gr4, + G4,, ) (2p + 1) 



( 4-36) 



which in turn implies that when vertical transport of angular momentum is taken into account, the mass accretion rate 
is enhanced by a factor (1 + G^z/Gr<p){2p +1). Conversely, wh en the vertical angular momentum flux is neglected. 

This result can be used to 



then T, 



r0 



;ed by a lactor [i + u^z/<~^r<p)[^P + i-)- t-^onverseiy, wn en tne vertical angular n 
scales as VrV^, a result obtained in numerical models i Hawlev fc Krolikll2002|) . 



estimate the minimum mass accretion rate at ri, where the MRI-driven turbulent stresses saturate at quasi-thermal 



levels (i.e. tr0(n) ^ p(ri)c^(ri)). Using Ma = 27rrl](- 



A/a(r) 



-Vr) gives 

n {-Vr, 



M, 



Edd 



( 4-37) 



where 

MEdd = 

e = O.le, 



= Io('^-r/mp)pdz 

AnGMuip/ieaTc) 



^SaT/nip is the Thomson optical depth of the disk over its vertical scaleheight and 
0.3eo iMyMgyr-i is the Eddington mass accretion rate for a conversion efficiency 
Q]^. The condition — 'Dr(n) — — ir0(n)/i^K('"i) — Cs(ri)/i)K(n) satisfied when angular momentum is solely 
transported by the turbulent r(f> stress then implies a dimensionless mass accretion rate 

rh,{rO^^^> 3 X 10-^ eo.iMrOrfilnin) , (4-38) 

MEdd 

where T5(ri) = T(ri)/10^K and T5(ri) = TT(ri)/10^ are approximately unity for a typical AGN disk. For ex- 
ample, accreting gas will come into thermal and radiative equilibrium and cool to the local blackbody tempera- 
ture Tbb ^ 2 X 10^(Ld/LEdd )^^'^(n/10?-s)~^'^^M7~^^'^ K when the number density reaches no{n) ~ IQi^ cm'^ (e.g. 
iKuncic. Celotti fc RBeslll997^ , which implies T^^in) ~ 10^ at n = lOr^ ~ 1.5 x lO^^^^ 

cm for a thin disk with 

h{ri) ~ 10~^ri. An immediate implication of this result is that since Afa is constant when there is no wind, then 
unless the vertical Thomson depths are unusually high, moderate to high mass accretion rates in geometrically-thin, 
optically-thick disks cannot be attributed to weak-field MHD turbulence alone. Indeed, substantially sub-Eddington 
mass accretion rates are a characteristic property of numerical models where angular momentum is transported solely 
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by tr<i, (see e.g. iHawlev fc BalbusI 120021 iHawlev et all 120011 iStone fc Pringldl200H) . This limitation arises from the 
weak-field nature of the MRI, which becomes ineffective once the turbulence reaches thermal levels. This constrains 
the maximum inflow speeds that can be attained in the accretion flow such that l^rl/cs — Cs/^k- Therefore, unless 
there is an additional instability capable of driving MHD turbulence from the weak to the strong regime, angular 
momentum transport in disks with moderate to high mass accretion rates must be largely attributed to a combination 
of vertical mass loss and mean-field torques (including large-scale MHD outflows, which we have not explicitly included 
here) . 

Note, however, that the role of the G^z torque in vertical angular momentum transport remains unclear. In princple, 
G^z could be larger than G^^ by as much as r/h, in which case this surface torque could play a key role in the energetics 
of thin disks. Consider the condition H4-31|) for angular momentum conservation. For a Keplerian disk, this implies 

M,-^ (tok) + ^ {2nr^Tr^) = ( 4-39) 

ar ar 

It is then straightforward to show that the (pz stress has an important effect on angular momentum transport when 

Z^III^ . 1 ( 4-40) 

and that the cjjz stress is important relative to the r(j) stress when 

r 



{tr4>)o ( 4-41) 



This extraordinary result reflects th e fact that {tct ,z)^ acts over the surface of the disk but (^^0)0 acts only over the 
height. The simulations reported bvl Miller fc Sto ne (2000) produce a volume-averaged {t^^)^ ~ 0.02 x (ir0)o- Even 
this low a stress is capable of physically interesting effects, as we argue in the following section when we consider the 
implications of mass and momentum transport for the energetics of turbulent accretion disks. 

5. THE DISK ENERGY BUDGET 

We now apply the above deductions from the momentum balance in a magnetized, turbulent accretion disk to 
calculate the power generated in the disk and emerging from the disk surface as both radiation and mechanical plus 
non-radiant electromagnetic energy. Dissipation of the latter two components may generate further emissivity in a 
corona and/or outflow. Our basis for this calculation is the total energy equation (|3-15() . Usually in expositions of 
accretion disk theory, the internal energy equation is used, since this explicitly identifies energy exchange processes, 
in particular, the transfer of free energy from the Keplerian shear in the bulk flow to the stresses which then dissipate 
that energy in the disk. In a turbulent MHD disk, this energy transfer is largely attributed to the turbulent energy 
production term {tij)sij ~ {tr4,)srij>, that is related to the rate of viscous dissipation via a turbulent cascade and the 
rate of bulk heating of the fluid via various additional radial and vertical transport terms. However, the internal 
energy equation for a turbulent, MHD accretion disk involves many intermediate terms that are cumbersome. The 
total energy equation, on the other hand, has the advantage that all terms are in total fluxes, which can be integrated 
straightforwardly to obtain a conservation equation in terms of the nett power in each of the available energy fluxes. 
We therefore prefer to use the total energy equation, which more elegantly describes the contributions of all terms to 
the overall energy balance in a disk. 

5.1. Order of magnitude estimates 

The total energy flux contains a number of terms and in estimating the relative importance of these, it is useful to 
deduce order of magnitude estimates of the various components of velocity. These are the mean velocity components, 
Vr, — {GM/ry/^ and v^, and the corresponding fluctuating components v',., v'^ and v'^. 

The component Vr is the radial inflow speed associated with the mass accretion rate viz. Vr — — Afa/27rrS]. An 
estimate of this quantity depends upon which stress dominates equation H4-39|) . If the rcj) stresses dominate, then 
\vrV^\ ~ |Tr0|/E, implying that 

\Vr\ _ Cs |(tr0)| ^ 
Cs Wk PoCs 

If the (jjz stresses dominate angular momentum transport, then approximating the radial derivative in equation (|4-31|l 
by division by r, we obtain 

\Vr\ r Cs 



Cs h Vk pqc': 



2 



5-2) 



If we take K^r^)! ~ Pov'^ (where v' is the magnitude of the turbulent velocity), then since the turbulence is weak 
{v'^ <Cg), the first estimate of Vr based on the rcj) stresses implies |?)r|/cs <C 1 for an optically-thick disk (see also the 
estimates in ii l4.5l above'). If the (t^z) stresses are of similar magnitude to (tr^) then Vr is larger by a factor of r/h when 
angular momentum transport by {tcf,z)'^ is important. Thus, if the disk is sufficiently thin that haM/r< (t^z)^ / {trtj>) , 
then the (jjz stresses could increase the radial Mach number Vr/cg to a value above Cs/wk- So far, simulations of 
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turbulent accretion disks fe.g. lMiller fc Stoii3l200CI() estimate these stresses to be approximately an order of magnitude 
smaller than the (tr^), since the MRI preferentially amplifies the r and 4> components of the fluctuating fields. 

The conclusion from these estimates is that the radial Mach number is likely to be small under most conditions. This 
means that most terms in the total energy equation involving Vr can be neglected, except those involving multiplication 
by a large quantity such as ~ or the gravitational potential. 

As far as the relative values of kinetic and magnetic stresses are concerned, we rely to some extent, though not 
exclusively on numerical results, which show that the turbulent magnetic stresses are a factor of a few higher than the 
Reynolds stresses. 

The fluctuating velocity components v'^ and v'^ are less than the sound speed, so that we neglect the turbulent energy 
flux terms related to these quantities. However, vertical gradients in turbulent energy flux terms involving v'^ should 
in general be retained since, in addition to being affected by the MRI, this component of the fluctuating velocity can 
also be enhanced by the Parker instability and the total rate of vertical energy transport can be nonnegligible in a 
thin disk. Turbulent perturbations in a magnetized disk become unstable to the gravitational modes of the Parker 
instability as matter drains down slightly elevated flux tubes, thereby enhancing their buoyant rise to the disk surface. 
The buoyant velocity field generated in this way can legitimately be considered as a turbulent velocity since regions 
of rising underdense flux tubes are balanced by falling, overdense regions and the nett mass-averaged velocity is zero. 
We estimate the buoyant velocity of such a flux tube by balancing the buoyant force per unit length corresponding 
to a density deficit of —6p with the drag force oc Cdv'J^ where Ca ~ 1 is the drag coefficient. Thus we estimate the 
vertical buoyant velocity of a tube of radius i?tubc at a height z to be (e.g. IParker 1979^1 : 



{My 



CdJ \ p J \ J \ r 



< 



X 1/2 /|c | \ 1/2 / n \ 1/2 



CdJ \ p J V 



CO ( 5-3) 



This velocity depends upon the density contrast generated, the height of the tube above the midplane and the tube 
radius, all of which are uncertain. For modest values oidp/p and i?tubc/^av the buoyant velocity can be an appreciable 
fraction of cq but it unlikely to be greater than cq. The Miller & Stone ( 2 000|) simulations, show an average kinetic 
energy associated with v'^ comparable to that associated with v'^. and v^, consistent with this. 

As we have indicated earlier, the vertical mean flow velocity at the base of the corona, w^, that could arise from an 
imbalance between the thermal and radiation pressures and the local gravitational field, cannot be estimated without 
a specific model. Nevertheless, the order of magnitude argument advanced in §0| suggests that it may be as large as 
the disk sound speed. Lacking any detailed knowledge of the magnitude of this velocity component, we retain terms 
in the vertical energy flux associated with it. 

In applying the total energy equation, we also neglect the viscous and resistive terms in the energy flux, — {iljv[j) -f 
As we noted in § 13 these terms could be locally important in, for instance, shocks or reconnection regions 
but when the energy equation is integrated over a large volume, their contribution to the surface integral is minor (see 
e.g. iHawlev & Balbus 2002, for a discussion on resistive heating in turbulent disks). Moreover, it is unlikely that the 
almost singular surfaces on which these terms are important would intersect substantially with the annular surfaces, 
r < r < r + dr, —h<z<h, over which we integrate the total energy equation. Also, neither term is important in the 
high wave number end of a turbulent cascade since they involve one less spatial differentiation than the Joule heating 
and viscous dissipation terms that appear in say, the internal energy equation. Those terms are locally important but 
they are not required in the integration of the total energy equation. 

5.2. Application of the total energy equation 

Let us now calculate the nett power available from each of the energy flux terms in the statistically-averaged equation 
for total energy conservation, given by H3-15|l . taking into account the order of magnitude estimates deduced in the 
preceding section. For an axisynimetric, steady-state disk: 

-^{r{F?)) + - , ( 5-4) 

r ar az 

where {F^) and {F^) are the radial and vertical components of the total mean energy flux (Ff), defined by H3-17(l . 
Integration over the disk height yields: 

^ j 2^r{F^) Az + 47rr(F,^)+ = ( 5-5) 

and this equation is the basis of the calculation of both the differential and total radiative energy flux from the disk. 
The dominant terms in the radial component of the total energy flux give 



1-2 



+ <t>G PVr - {tr4,)v4, + {Fr) + {Qr) ■ ( 5-6) 
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Upon vertical integration over a thin disk, followed by differentiation with respect to radius the last two terms become 
0{ha,v/r) compared to analogous terms in the vertical energy flux, implying 



1^2 



where T^^ is the vertically- integrated r(j} stress defined by (|4-30|l . The term —2T:rTr^v^ 
using equation H4-33() for angular momentum conservation: 



(5-7) 

2'KrTrtbVK can be eliminated 



dr 



dr' 



where 



1/2 



The vertical component of the energy flux is 
'1 



E\+- 



+ {ph + {u^v^y + (mbWz) 



,/\ + 



(5-8) 
(5-9) 

( 5-10) 



There are some terms in this expression such as the first two terms representing the vertical flux of kinetic plus 
gravitational energy in a wind and the term (t^z)^^^, that we should obviously keep. There are also some terms in 

this equation that we can immediately discard. The term p+/i+'D+ = (m+ + p'^)v'^ is small at the photosphere-corona 
boundary. For similar reasons we can also neglect the turbulent enthalpy flux term {ph'v'^). In view of the order of 
magnitude estimates considered in § 15.11 the remaining terms may be comparable. To simplify the final expression, 
we combine the remaining terms involving a vertical velocity component into an energy flux term {'>pz)~^ , defined by 

With these simplifications, the vertical component of the total energy flux reduces to 

^ ' ^ ^ ■ '■ + 47rr((V'.)+ + (Fz)+ + (Qz)+) . (5-12) 



( 5-11) 



47rr(F. 



E\ + 



47rr —V, 



p+vt 



Subsitution of equations H5-7|l and H5-12(l into the disk total energy equation H5-5|l . utlizing equation H5-8|l gives the 
following for the radiative flux from both sides of the disk (after some algebra): 



3GMM 



2 r2 
2 



dr' 



r j^. [2 

-47rr(V'.)+-47rr(gz)+ . (5-13) 

The first term on the right of this expression represent s the binding energy fiux associated with the nett mass fiux and 
is familiar from standard disk theory (c.f.lShakura & S unvaevll973() . The second term represents a correction resulting 
from a non- vanishing stress at the inner radius. The third term is related to the flux of kinetic and gravitational energy 
in the wind; this term is discussed further below. The fourth term represents the effect of the (jfz turbulent stresses 
and its magnetic component is intimately associated with the Poynting flux into the corona; this term is also discussed 
further below. The fifth term, — (V'z)"'', represents another component of the wind power associated with the advective 
part of the Poynting flux and advection of other turbulent quantities. The fifth term, —{Qz), represents the effect of 
an external heat flux arising from above the disk. Since (Qz) < 0, this heating enhances the radiative flux from the 
disk. 

The above form of the radiative flux is convenient for integration since it involves the (constant) total mass accretion 

rate, M in the first term and the factor CK^^(r) has a simple analytical form. Nevertheless, in order to facilitate 
comparison of the various terms, the third term may be integrated by parts and partially combined with the first term 
to give: 



+ 3GMMa (2) 

2 



2 r 



3TTr^Tr.^{n)- 

r 



where 



d'^(0 = 1 



n^i/^ A/a(n) 

Ma(r) 



( 5-14) 
( 5-15) 
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The form H5-14|) of the radiative flux from the disk more clearly shows the effect of the local accretion rate in the first 
term and the effect of the angular momentum transported by the wind and the (pz stress in the third and fourth terms 
respectively. 

An expression for the total disk radiative luminosity, id, may be obtained by integrating each term of equation (|5-13|l 
over all disk radii, from r — io r — oo. This also helps to elucidate the significance of the various terms in the 
expressions for the radiative power from the disk. The result is: 

IGMM ^_ 2 /„^ IGAfMw(ri) /"^ 1 GA/ dM« 



id =2^— - 27rr,T.,(n)f^K(n) - |^ 

A'Kr{t4,^{r)) + v^{r)dr - I 47rr (-02)+ dr - / 47rr (Q)+ dr (5-16) 



The first and third terms combine to give the familiar expression for the accretion power: 

2 n 

In conventional accretion disks, this is the dominant term responsible for the disk luminosity. In the generalized 
formalism developed here, however, we take into account disk solutions in which Ma decreases towards small r owing 
to the vertical loss of mass. Thus, for disks with strong wind mass loss, Ma(n) is a fraction of the nett mass flux M 
and Pa is not necessarily the dominant term. 
In considering the effect of a nonvanishing stress at rj, we introduce the parameter 

Outward angular momentum flux due to turbulent stresses at r = r\ 



Inward angular momentum flux due to accretion at r = ri 



( 5-18) 



Note that since Tr^{ri) < 0, consistent with the stresses defined elsewhere in the disk, this parameter is positive. 

We now combine the fourth, fifth, and sixth terms on the right hand side of H5-16|l into a single term describing the 
total power removed from the disk by a wind: 



Pw = 



1 GM 
2~V 



47rrp+w+dr- / 4:nr(t^^{r)) + VK{r) dr + / Airr {^p^)+ dr . (5-19) 



The leading term in the wind power is negative since the disk material is highly bound; the only way in which the 
wind power can be positive is for the succeeding terms in to be positive. (See the discussion concerning winds 
below equation 1)4-211) '). 

We define the nett heating rate of the disk due to irradiation and heat conduction from an external corona by 

/•oo 

Fq = -/ 47rr(Q)+dr>0 (5-20) 



We can now write the total disk luminosity as the following sum: 

id - (l + 2K)Pa-Pw+PQ . (5-21) 

When vertical energy transport from the disk (i.e. P^) is important, the disk luminosity is reduced as a result of the 
conservation of energy. Realistically, some of the total power P^ removed from the disk is available to heat the corona 
and potentially, the most important contribution to this term is: 



4^r-(t^,)+WKdr , (5-22) 

which represents the integral of part of the Poynting flux over both surfaces of the disk (see equation 1)3- 18() ). To 
evaluate the importance of this term we compare the integrand of 1)5-22)) to the leading term in equation ()5-14)) for the 
disk energy flux. This comparison shows that P^^ is an important component of the disk energy budget when 



1 ( 5-23) 



3 47rAfaWK 

We showed in §0](see equation (I4-4UII ') that (t^z) has an appreciable effect on the accretion process when 

-J^tMl . 1 . ( 5-24) 

Mal^K 

Hence, when this is the case, the power P^^ is comparable to Pa, that is, when the 02; stresses are important for 
angular momentum transport, they do significant work on the flow. 
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Of course, the effect of the cjjz stress would be unimportant if it fluctuated in sign. However, if the disk has a 
wind, then one expects {t(f,z)^ < as a result of inertia causing magnetic field loops to trail the rotation of the 
disk. When (t^z)"*" is systematically negative, the disk luminosity is reduced for a given mass accretion rate and a 
corresponding amount of energy is available to heat the corona. As we have shown in § 01 even a weak stress can be 
dynamically important and the above shows that it can also be energetically important. I ndeed, the importance of (jj z 
stresses in thin disks is widely appreciated in accretion-outflow models. In the model bv lPudritz fc NormanI ijlDS,*^ . 
for instance, the entire accretion flow is attributed to outward angular momentum transport by the torque associated 
with these stresses. A major physical difference between our model and that of disk-driven wind models such as that 
of Pudritz & Norman ( 1983) and Blandford & Pavnc ( 1982) is in the mean magnetic flux threading the disk. 

We can also express the condition (|5-23|l in terms of other disk parameters using. A/a ~ 47rr/iavPo|'5r |- For P^z to be 
a significant fraction of the accretion power 



PO /lav Vr Po 



For representative AGN values po / p^ 



-i6 



0+ 



r Vk p^ \ r J Co 
10'', hav/r ^ 10^'', this condition effectively implies (c.f. eqn. 

•"l B' Co 



( 5-25) 



( 5-26) 



If the radial inflow velocity is subsonic, and there is at least a modest vertical field (e.g. B'^/B'^ ~ 0.1, say), then 
this condition may be satisfled if the coronal Alfven speed is comparable to the Keplerian speed. However, depending 
upon the details of disk parameters (including h^-^ jr and the radial Mach number), the condition may be satisfled for 
much lower values of the Alfven speed. 

At this stage we also consider further the relationship of the power P^z, to the vertical component of the Poynting 
flux: 

(^z) = («b)5z + (j^bO - {tlXv, - if^Xv^ - - (t^?z<> - (iL^^;> - ^f.A) ■ ( 5-27) 



The term 



( 5-28) 



in this expression contributes directly to P^z. (Other terms in (S'z) are represented in -i/iz-) The power P^z represents 
the integral of this part of the Poynting flux. It is interesting that the term (t^^)^^ does not appear in the expression 
for the local radiative flux, viz. equation H5-13|l . Such a term is initially present as the above derivation for the local 
radiative flux shows. However, it is eliminated by a corresponding opposite term arising from the expression for T^^. 
What does remain in the expression for the local radiative flux is a more complicated term involving multiplied by 
an integral over (t^^). This cause of this is that the angular momentum flow is driven by the total stress whereas the 
Poynting flux is only related to the magnetic part of that stress as well a the inter-realtionship between the rt/) stress 
and the (^z stress in the angular momentum equation 14-31|) . 

The term Airripz dr involves the integration of a number of terms all of which could be important in the disk 
energy budget. Rather than consider them all here, we consider the representative, and probably most important 
terms: 



47rr(5'f)) =47rr((wB - ^zz)^'z)"^ = 47rr 



{B'/ + B'/)+^ 

'a ^ 

47r 



( 5-29) 



( 5-30) 



(2) (3) 

The term {Sz ) represents the contribution to the Poynting flux from the wind; the term {Sz ) 
contribution to the Poynting flux from turbulent diffusion. Consider the flrst term. 



4:TTp + 



Comparing this with the leading term in equation 

(vl) + \dAU/dr\ (vj 



dM^ {B'/ + B'/) + 
dr 47rp+ 
2} for 47rr(Pz) we have: 



dM^ 
dr 



represents the 



( 5-31) 



3GMMj2r^ 



r 



dM^ 



dr 



r dMa 
Ma dr 



( 5-32) 



For the power associated with this term to be important, it would seem that both the wind mass loss rate should be 
comparable to the mass accretion rate (r|dAfw/dr| ^ Afa) and the Alfven speed in the corona should be comparable 
to the local Keplerian speed. Otherwise, the requirements on each of these factors would probably be excessive. 
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Similarly, one can compare the diffusive terms to the same leading term in AiTr{Fz). The diffusive flux is important 
when 

.(3)\ 

1 ( 5-33) 



iGMM/2r^ I GMM^ 

This condition is similar in some respects to the one above referring to the systematic velocity component. However, 
in this case the turbulent velocity is not related to the wind mass-loss rate and one should allow for the possibility 
that the turbulent velocity can in principle be larger than the mean wind velocity. Again, utilizing Ma ~ 47rr/ipo|^r| 
in this expression leads to the condition 

— 2—-^ r (5-34) 

As in the treatment of the P^z term above, typical AGN values po/p^ ^ 10^, hav/r ~ 10"'^ this condition relies on the 
ratio of the radial inflow velocity to the turbulent buoyant velocity. Nevertheless, the impression from (|5-34|) is that 
the Alfven speed in the corona should be fairly high with respect to the local Keplerian speed for turbulent diffusion 
to be important. 

5.3. Relative importance of Poynting flux components 

We have seen above that the importance of the various terms in the Poynting flux compared to the accretion power 
generally depends upon the ratio of the coronal Alfven speed to the local Keplerian speed and its relationship to disk 
and wind parameters. The above relationships give us a good idea of the importance of the various components of 
the Poynting flux in absolute terms. Further insight is obtained if we compare the various terms with respect to one 
another; this informs us of the conditions under which each component is likely to dominate. In these comparisons 
we remind the reader that 5*1^^ is associated with the power P^z, si^"^ is associated with wind-advective part of the 

Poynting flux and Sz is associated with the turbulent diffusive part. 
The ratio 

( 5-35) 



-{B'^B'^)+Vk -{B'^B',)+ r Co 

Since Vz ^cq, the vertical magnetic field would have to satisfy B'^/B' <hay/r <C 1 in order that the wind advective 
power dominate Ptj,z- 

Similarly, the ratio of the turbulent diffusive power to the P^^ is determined by the ratio 

^(3) Ib:1w\+ /n'2\ „ ,, 

( 5-36) 



S'z'^ (P'2)+ K 



Si'^ -{^)+v^ -(B',B'^)^ r CO 

Again, given that v'^ < cq the vertical component of magnetic field would have to be extremely small for the turbulent 
diffusive term to dominate. 

5.4. The limiting case of vertical transport only 

In a real disk, angular momentum is probably transported by a combination of both radial and vertical fluxes. The 
limiting case of radial transport only has, of course, been well-explored in the past. Let us now investigate the opposite 
limiting case of vertical transport only, specifically as it affects the estimate of the power of the disk wind. In this 
limiting case, equation (|4-31|) becomes: 



_d 

dr 

In this limit, the non-kinetic, non-gravitational part of the wind power is dominated by P^z so that 



M.v^r = v^r^ - ^nr^t^z)^ ( 5-37) 



p 



2 r dr 

Using equation H5-37|) the Poynting flux term can be expressed as 



dr ( 5-38) 



- 4^r(i^,)+«K = (™k) ( 5-39) 

r dr 

Inserting this expression into the integral H5-38|) for the wind power, one obtains: 

IGMM^ (5-40) 

2 ri 

i.e. the wind power is equal to the accretion power. The corresponding disk luminosity is 

id ^ Pq ( 5-41) 
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i.e. the entire disk luminosity is attributed to external feedback heating from the corona. In this limiting case, some 
of the wind power would be dissipated in the corona and the rest would escape to infinity. 

There is a difference between the power in such a wind and the pow er in a centrifug ally driven wind with a nett 
magnetic flux threading the disc. For example the power in afBlandford fc Pavnd l)1982|) wind is 

Pw.BP oc Bl{n)rl f ^ j ( 5-42) 

and depends upon the nett magnetic flux in the disk, as well as the Keplerian velocity at the inner radius. 

In the theory we have developed here, the nett flux is zero, so that the magnetic field plays a different role than that 
envisaged in centrifugally driven flows. Since the r(j) stresses play no role, the result represented in equation (|5-40() is 
inevitable. The binding energy released by the accreting disk material has to be manifest in the power of the wind. 

6. SUMMARY AND CONCLUSIONS 

In this paper, we have established a self-consistent framework for the theory of magnetized, turbulent disk accretion 
around black holes. Our formalism is the first to consistently treat turbulent disks using a robust statistical averaging 
procedure that explicitly includes dynamical equations for the evolution of the magnetic field together with the con- 
servation equations for mass, momentum and energy transport. We have paid special attention to vertical transport 
of conserved quantities and consistently related such transport to the dynamical structure of the underlying disk. 
Although the nett magnetic flux is assumed to be zero, the formalism is nonetheless sufficiently general to allow the 
straightforward inclusion of a systematic nett mean- field component. 

Our main results can be summarized as follows: 

1. We have derived a comprehensive set of equations describing the transport of mass, momentum, internal energy, 
turbulent kinetic and magnetic energies, and total energy. The statistically- averaged conservation equations are 
completely general and are applicable to other subject areas involving turbulent magnetic fields. 

2. We have applied the statistically averaged equations to a geometrically-thin, optically-thick accretion disk that 
is stationary and axisymmetric in the mean. We have demonstrated that when vertical transport of mass, radial, 
vertical and angular momentum, and energy is self-consistently treated, the general equations include additional 
terms related to a disk wind and turbulent azimuthal-vertical stresses on the disk surface. 

3. We have shown that the total azimuthal-vertical stress can have a significant dynamical and energetic effect on 
the disk even though it may be numerically small compared to the radial-azimuthal stress that has dominated a 
large amount of accretion disk theory and simulation, t o date. Note, however, t hat the importance of this stress 
has also been realized in accretion-outlfow models fsee lKonigl fc Pudritdl2?KKi for a review). 

4. We have derived an expression for the radiative luminosity from the disk photosphere and shown clearly how 
this relates to the mechanical power in a wind and to the Poynting flux, thereby identifying the possible sources 
responsible for powering coronae and/or outflows from accreting black holes. This expression also entails a 
different distribution of radiative flux than a standard accretion disk. This in turn affects the integrated spectrum. 
Again, we defer the details to future work. 

5. We have discussed the three main sources of Poynting flux into the corona - a component associated with the 
product of the azimuthal - vertical component of the turbulent magnetic stress, a component associated with 
wind advection of magnetic energy and a component associated with turbulent diffusion of magnetic field from 
the disk into the corona . The first component probably dominates in most cases even if the azimuthal - vertical 
stress is quite small in comparison to the radial - azimuthal stress. In the course of the analysis of the condition 
for a wind and the conditions for a significant Poynting flux into the corona, the ratio of the coronal Alfven 
emerges as a critical parameter. When this ratio is of order unity, important magnetic effects are clearly present. 
However, there is also the prospect of significant effects when this parameter is less than unity. This region of 
parameter space is currently a relatively unexplored avenue of research in black hole accretion disks. 

6. In the limiting case, when all of the angular momentum transport is through the vertical-azimuthal stress, we 
have shown that the wind power, at the base of the wind, is exactly equal to the accretion power. Some of 
this power would be dissipated in the corona. This is the first time that a coupled disk-corona model has 
demonstrated, in a physically consistent fashion, the possibility of significant power emanating from the corona. 

7. This is also the first time that the power of a disk wind has been dynamically linked to the process of accretion. 
In models with a nett magnetic flux the wind power is related to the strength of the magnetic field. 

8. The existence of a coronal wind and the production of intense coronal emission are inextricably linked. The major 
influence in heating the corona is the stress that is responsible for transporting angular momentum vertically. 
The wind is essential to transport this angular momentum away from the disk. 
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In a realistic disk, we expect that both radial-azimuthal and aziniuthal-vertical stresses would be involved in 
the transport of angular momentum, as well as a nett mass loss from the innermost regions. Nevertheless, our 
limiting solution provides a good physical basis for the commonly held notion that accretion power could be channelled 
into significant coronal emissivity and/or outflow comprised of electromagnetic and bulk kinetic components. Thus 
there is the prospect of explaining not only the coronal emission from radio-quiet AGN and Galactic Balck Hole 
Candidates (GBHCs), but also systems such as Ultra-Luminous X-ray (ULX) sources and analogous AGN sources 
such as Broad Absorption Lin e (BAL) quasars where large mass outflows (e.g. M^dd) are inferred (see for 

example iKing fc Poundsll2003l and references therein). 

With a self-consistent framework now established, we are in the position of being able to consider further, via 
specific models, the complex relationship between disk, corona and outflows in a variety of sources, and to examine 
more thoroughly the conditions for the initiation of a wind and the implications for the general structure of the 
immediate environment of accreting black holes. 
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